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Abstract 

An acoustic oscillation of the primeval photon-baryon fluid around the decoupling 
time imprints a characteristic scale in the galaxy distribution today, known as the 
baryon acoustic oscillation (BAO) scale. Several on-going and/or future galaxy sur- 
veys aim at detecting and precisely determining the BAO scale so as to trace the 
expansion history of the universe. We consider nonlinear and redshift-space distor- 
tion effects on the shifts of the BAO scale in /c-space using perturbation theory. The 
resulting shifts are indeed sensitive to different choices of the definition of the BAO 
scale, which needs to be kept in mind in the data analysis. We present a toy model to 
explain the physical behavior of the shifts. We find that the BAO scale defined as in 
Percival et al. (2007) indeed shows very small shifts (^ 1%) relative to the prediction 
in linear theory in real space. The shifts can be predicted accurately for scales where 
the perturbation theory is reliable. 

Key words: cosmology: large-scale structure of universe — theory — methods: 
statistical 



1. Introduction 

The baryon acoustic oscillation (BAO) is an oscillation of photon-baryon fiuid imprinted 
in the matter spectrum as a characteristic signature. Recently it was detected in the SDSS 
and 2dFGRS galaxy distribution (e.g., Eisenstein et al. 2005; Hiitsi 2006; Tegmark et al. 2006; 
Padmanabhan et al. 2006; Percival et al. 2007; Cole et al. 2005), while its counterpart in 



the cosmic microwave background (CMB) has already played an important role in precision 
cosmology (Spergel et al. 2007). Currently, using the BAO scale as a standard ruler is regarded 
as one of the most promising tools to trace the cosmic expansion history. This characteristic 
scale basically corresponds to the sound horizon at recombination (e.g., Eisenstein et al. 2005, 
and see also Appendix A): 



where Cs{z) is the sound speed at redshift z, and Zj-^c is the redshift at recombination 1089). 
In the second equality, keq is the horizon scale at the matter-radiation equality epoch, Zeq, 
i?rcc = R{^rcc) a^d R^q = R{zcq) are the ratio of the baryon to photon momentum densities 
at z^ec and Zf,q. Finally the last equality is an approximate fit where Qm and Qf, are the 
density parameters of matter and baryon, and h is the current Hubble constant in units of 
lOOkms-^Mpc-^ 

The BAO length scale itself can be computed accurately and is indeed insensitive to the 
presence of dark energy that affects the expansion of the universe at relatively low redshifts. 
Precisely for these reasons, it is a useful standard ruler of the universe. In particular, it is 
supposed to be a good tracer of wbe =Pde/pde, where pbe and pde are the pressure and the 
density of dark energy (e.g., Blake & Glazebrook 2003; Seo & Eisenstein 2003). Also it can 
potentially be used to falsify the possibility of alternative law of gravity which might explain 
the acceleration of the cosmic expansion (Shirata et al. 2005; Yamamoto et al. 2006). 

Consider a sample of galaxies with measured redshifts, the observed BAO scale provides 
estimates of the angular diameter distance Da{z) and the inverse of the Hubble parameter 
1/H{z), which correspond to the scales perpendicular and parallel to the line-of-sight direction, 
respectively. They in turn can be translated into the estimate of wbe- Figure 1 shows how the 
fractional errors of three important scales, the angular diameter distance Da{z), the inverse of 
Hubble parameter 1/H{z) and their average over three dimensions {D\{z) / H{z)Y^^ propagate 
to that of wbe- The cosmological parameters assumed are the third year WMAP results (Spergel 
et al. 2007, see section 3 for details). The two shaded regions show the approximate targeted 
redshift ranges of a future galaxy redshift survey, WFMOS (Wide-field Fiber-fed Multi-Object 
Spectrograph). Typically a ratio of Aw/w and Ad/d around z = 1 ranges from 3 to 5, while 
the value depends slightly on a specific choice of cosmological parameters. Thus the ~ 3% 
determination of wbe requires the sub-percent accuracy/precision in determining the BAO 
scale, which is challenging from observational, and even theoretical, points of view. 

Until recently it was often assumed that gravitational non-linearity and redshift- 
distortion effects do not shift the BAO scale, although they significantly affect the amplitude of 
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Fig. 1. The error propagation from measured scales, d, to the dark energy equation of state 
parameter, wde, as a function of redshift. We choose 1/H{z) (dotted), and Da{z) (dashed) 
for d, which correspond to the separations parallel and perpendicular to the line-of-sight di- 
rection. We also plot the three dimensional average, {D\{z)/H{z)y/^ (solid) for d. The 
shaded regions indicate the targeted redshift ranges of a future galaxy survey, WFMOS. 

the BAO. Some simulations (e.g., Meiksin, White & Peacock 1999; Seo & Eisenstein 2005; 
Springel et al. 2005; Angulo et al. 2005; Jeong & Komatsu 2006; Ma 2006; Eisenstein 
et al. 2006; Angulo et al. 2007) and analytical works (e.g., Eisenstein, Seo & White 2006; Guzik, 
Bernstein & Smith 2007; Smith, Scoccimarro & Sheth 2007a; Smith, Scoccimarro & Sheth 
2007b) have been done to study these effects on the power spectrum or two-point correlation 
function at the BAO scale. For a few percent accuracy on u'de? more precise BAO predictions 
are required, which we consider in the present paper on the basis of the one-loop corrections 
from perturbation theory. 

Very recently, several authors (Crocce & Scoccimarro 2006a; Crocce & Scoccimarro 
2006b; McDonald 2007; Valageas 2007; Padmanabhan & Ray 2006; Matarrese & Pietroni 2007a; 
Matarrese & Pietroni 2007b; Crocce & Scoccimarro 2007) attempt to take account of higher- 
order corrections. For the first time, we take account of the redshift-space distortion effects on 
the BAO scale applying the result of Scoccimarro (2004). 

The outline of this paper is as follows; section 2 briefly summarizes the perturbation 
theory of density and velocity fields. The details of the numerical calculation are shown in 
section 3. The results are shown in section 4 with discussion and interpretation in section 5. 
Finally, section 6 is devoted to the conclusion of this paper. 



2. Non-linear power spectrum in real and redshift spaces 



Cosmological perturbation theory predicts tlie gravitational clustering of matter distribu- 
tion by a systematic expansion of the cosmological density and velocity fields. It provides us an 
accurate prediction to the non-linear modification of matter power spectrum as long as the grav- 
itational clustering is mildly non-linear (e.g., Juszkiewicz 1981; Vishniac 1983; Fry 1984; Goroff 
et al. 1986; Suto & Sasaki 1991; Makino, Sasaki & Suto 1992; Jain & Bertschinger 1994; 
Bernardeau 1994; Matsubara 1995; Scoccimarro et al. 1998; Chodorowski & Ciecielag 2002, 
and see Bernardeau et al. 2002 for a review). In this section, we briefly review the formalism to 
calculate the non-linear correction to the matter power spectrum. A model of redshift-space dis- 
tortion is also presented, which will be used later to estimate the non-linearity of redshift-space 
power spectrum. 

2.1. Perturbation theory 

To deal with the gravitational clustering of the matter distribution, we adopt the hy- 
drodynamic description and treat the dark matter and baryons as a pressureless dust fluid. 
According to Scoccimarro (2001), the evolution equations for the cosmic fluid can be written 
in a compact form by introducing the two-component vector: 

-S))- 

where the subscript i = 1,2 selects the density or velocity components, with 5(k) and 0{'X]ri) = 
V-v(x;?7) respectively being the Fourier transform of the density fluctuation and peculiar 
velocity divergence. The variable t] represents the time variable defined hj rj = lnD+, where 
is the linear growth factor normalized to unity at present. The quantity Ti. is the conformal 
expansion rate given by 7i = dlna/dr, where a is the cosmic scale factor, and the quantity / is 
the logarithmic derivative of the linear growth factor, /(?]) = d\nD^/d\na. 

Assuming the irrotational fluid flow, the evolution equations can then be written as 

dr,^i{k;r])+nij{r])^j{k;7]) = J d^hd^k2%jkiK'i^iM)'^ji^i;v)'^kik2;7]), (3) 

where the kernel 'jijk is the vertex matrix which represents the non-linear interaction between 
different Fourier modes: 
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with 6d being Dirac's delta function. The matrix Qij is given by 
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This expression is not exact in the non Einstein-de Sitter cases, but it still provides an accurate 
prescription. 

The evolution equation given above is systematically solved by a perturbative expansion 
of two- component vector. Ignoring the decaying modes, one obtains 

oo 

vi/,(k;r/) = ^e"''#)(k) 



(6) 



n=l 



(n) 



l2 



l3n— 3 
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k„)^(")(ki,---,k„)5i(ki)---5i(k„), (7) 



where the quantity 5i stands for the initial density fluctuation, which we assume is a Gaussian 
random variable. The Fourier kernel jFf"^(ki, • ■ ■ , k„) represents the mode coupling between 
different Fourier modes originating from the non-linear interactions. In Appendix B, the explicit 
expressions for J^^"^ are presented up to the third-order of perturbative expansion. 

We are interested in the non-linear evolution of power spectrum, which are evaluated as 
the ensemble average of the quantity 

(*,(k;r/)vl>,(k';r/)) = (27r)35^(k + kOP.,(A;;r/) = 1,2). (8) 

Note that we obtain the three different power spectra: P^^ from («, j) = (1, 1), P^q from (z, j) = 
(1,2) and (2,1), and Pee from (z,j) = (2,2). Substituting the solutions (6) up to the third-order 
perturbations into the above, the next-to-leading order corrections are obtained, which can be 
summarized as 



P,,{k-T^) = DlP'^{k) + D 



(9) 



In the above expression, the first term in the right-hand-side is the linear power spectrum given 
by 

(5i(k)(5i(k')) = {2nf8r>{\^ + \^)P\k). (10) 

On the other hand, the terms in the bracket of equation (9) are the so-called one-loop corrections 
as a result of the non-linear mode-coupling: 



p(13) 
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J^i'^ (k, q, -q) + (k, q, -q)] P'^{k)P\q) , 
'J^t\k - q,q)^f (k - q,q)l pL(|k - q|)P^(g). 



(11) 
(12) 



Note that the part of the integrals over the Fourier mode q are analytic, which we use. The 

(22) 

resultant expressions include the two-dimensional integrals over r and x for the P-^j {k) part, 
while the P/^- {k) part has the one- dimensional integral over r. Their explicit functional forms 
are summarized in Appendix B, together with the definition of r and x. 
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2.2. A Model of redshift- space distortion 

The observed distribution of galaxies constructed from redshift surveys is inevitably 
distorted due to the peculiar velocity of each galaxy. This effect, known as redshift-space 
distortion, is classified in two ways. On large scales, the bulk motion falling into a cluster 
apparently squashes the matter distribution, which enhances the clustering signal along the line- 
of-sight. On small scales, on the other hand, the virialized random motion of galaxies residing 
at a cluster suppresses the amphtude of the clustering signal. This is called the fingers-of-God 
(FOG) effect. 

Based on the linear perturbation theory. Kaiser (1987) proposed a formula for redshift- 
space power spectrum on large-scales: 

P^^\k,fi;v)=[l + fiv)f^fP^'\k;v), (13) 

where fi is the cosine of the angle between the line-of-sight direction and the Fourier mode 
k. The spectra P'^^\k]ri) and P^^\k, ij.;r]) respectively denote the matter power spectra in real 
and redshift spaces. Several authors proposed models for redshift-space power spectrum taking 
account of the small-scale random motion (Peacock & Dodds 1994; Park et al. 1994; Cole, 
Fisher & Weinberg 1994; Balhnger, Peacock & Heavens 1996; Magira, Jing & Suto 2000). In 
their models, the FOG effect is expressed by a damping factor, DpoG^k,^), which they assumed 
Gaussian or Lorentzian. Their models are written as 

P^^\k,^^;v)=[l + fiv)^^fP^'\k■,v)DFOG{k,^I). (14) 

Note that the damping factor DpoQ{k,fi) asymptotically approaches unity in the large-scale 
(small k) limit. 

We should notice that the models given in equation (14) have been constructed to deal 
with a relatively small-scale clustering. It might not be accurate enough for our interest in 
the precision measurement on BAO scales. As for the models (13), it has been advocated that 
nonlinear random motion can not be negligible even in the large-scale limit and one could not 
recover the Kaiser formula (13) (Scoccimarro 2004). We therefore look for alternative models 
relevant for the BAO scales. 

Scoccimarro (2004) recently proposed a physically plausible model of redshift-space dis- 
tortion based on perturbation theory. He improved the Kaiser formula to take account of the 
non-linear evolution of density and the velocity fields, as well as the FOG effect. Using the three 
different power spectra defined in equation (9), the explicit expression for the redshift-space 
spectrum becomes 

p(^)(A;,/i;r/) = [P5,(A;;r/) + 2/(r/)/i2p5,(A;;r/) + /(r/)Vne(A;;r/)^ 

xexp[-/(r/)V^V„2(^)], (15) 
with the quantity cr^(?7) being the one-dimensional hnear velocity dispersion: 
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The model (15) properly accounts for the non-linear mode- couplings of density-density, density- 
velocity and velocity-velocity fields in the Kaiser formula. In this paper, we adopt the model 
in equation (15) to calculate the redshift-space power spectrum. Note that our current model 
of the FOG effect is still empirical and has to be tested with numerical simulations, which will 
be discussed in future work. 

To compute P^^\ we use the one-loop results of power spectrum except for the linear 
velocity dispersion (16). The results are then presented by taking the angular average: 

P^\k;v) = l£^df,P^''\k,f,-rj). (17) 
3. Details of calculation 

3.1. Initial condition 

In what follows, we use the CAMB code (Lewis, Challinor & Lasenby 2000) to calculate 
the transfer function for the linear power spectrum P^{k). We assume a flat ACDM model 
with adiabatic Gaussian fluctuation. We use the best-flt values of the cosmological parameters 
determined from the three year WMAP data (Spergel et al. 2007): ^mh^ = 0.1277, ^bh^ = 
0.02229, h = 0.732, = 0.958, r = 0.089 and cig = 0.761, where is the scalar spectral index 
(without running), r is the optical depth and as is the rms of the density contrast smoothed 
with an 8h~^Mpc top-hat window. Note that from equation (1), we have r<, = 148Mpc, which 
leads to the oscillatory behavior like sin{krs) in the transfer function of the matter fluctuation. 

3.2. Accuracy of the perturbation prediction 

In general, the prediction from perturbation theory eventually breaks down when the 
non-linear correction dominates over the linear theory prediction. While we cannot rigorously 
deflne the validity range of the perturbation prediction, Jeong & Komatsu (2006) recently 
showed that the one- loop correction to the power spectrum Pssik) accurately describes the 
N-body simulations to better than 1% accuracy when 



A\k-z)^^^^^<0.4. (18) 

Thus, for the limitation of the perturbation predictions, we deflne the maximum wavenumber 
through 

A\k,%{z)-z) = 0.4. (19) 

Note that the limitation for the power spectra Pse and Peg would not be described by equa- 
tion (18) (see Crocce & Scoccimarro 2006b). Nevertheless, for simplicity, we adopt the above 
criterion in both real and redshift spaces. 

Figure 2 shows the numerical values of ki% as a function of redshift. Here, in addition 
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to the result for our fiducial cosmological model (solid), we also plot the results for the cases 
with slightly different amplitudes: ag = 0.7 (dashed) and 0.9 (dotted), keeping the other re- 
maining parameters fixed. The validity range of perturbation theory has a strong dependence 
on the redshift and the normalization. For the higher fiuctuation amplitude with = 0.9, the 
perturbation theory breaks down at a relatively smaller wave number. 




0.1 0.2 0.3 0.4 
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Fig. 2. The maximum wavenumbcr of the vaHdity range for the perturbation theory, ki%, defined by 
equation (19). The sohd fine represents the result for our fiducial model. We also plot the results for 
the cases with a slightly larger amplitude (erg = 0.9; dotted), and a smaller amplitude {as = 0.7; dashed). 



3.3. Characterizing the acoustic oscillation scales 

We are especially concerned with the systematic infiuences of the non-linear clustering 
and the redshift-space distortion on the characteristic scale of BAO, affecting its use as a cosmic 
standard ruler. To quantify the degree of these influences, we must flrst deflne the algorithm 
to characterize the oscillatory features from the power spectrum, by which the locations of 
peaks and the troughs are identifled unambiguously. Then, the systematic influences can be 
investigated by measuring both the amplitude and the locations of peaks (and troughs). There 
are several approaches to separate the oscillatory pattern from the power-law behavior of power 
spectrum. In this paper, we study the following three methods: 

(i) Divide the matter power spectrum by a smooth linear spectrum: 

fiUk) = ^y (20) 

where Pnw{k) represents the "no- wiggles" approximation of the linear power spectrum 
given by Eisenstein & Hu (1998), which can be evaluated for our flducial cosmological 
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model. With this characterization, smooth power-law behavior in the power spectrum is 
effectively eliminated and the peaks and the troughs of acoustic oscillations are clearly 
identified. Note, however, that the above definition only corrects the power-law feature 
of the linear spectrum and non-linear corrections of power spectrum are left untouched in 
the numerator. Furthermore, there is uncertainty in the normalization of the amplitude. 

(ii) Take the logarithmic derivative: 

/So«^^. (21) 

In contrast to the method (i), this method is free from the uncertainty in the normalization 
of amplitude. Also, it automatically separate the smooth power-law feature of P{k) from 
the oscillatory behavior. In practice, however, it seems rather difficult in estimating the 
logarithmic derivative from the noisy binned data. 

(iii) Divide the matter power spectrum by the smooth spectrum constructed from the matter 
power spectrum itself: 

= (22) 

smooth l^"' J 

This method is almost identical to the one proposed by Percival et al. (2007). Since we 
do not need a reference spectrum, it seems useful in characterizing the acoustic signature 
from the observed spectrum. Here, the smooth spectrum is constructed as follows. First 
we sample 70 data points from our perturbation spectrum between k = 0.02/iMpc~^ and 
0.3/zMpc~^ with equal separation and without weight. Then we select 8 nodes separated 
by Ak = 0.05/iMpc~^ between 0.025/iMpc~^ and 0.375/iMpc"^, and an additional node 
ai k = 0.001/iMpc"^. Connecting the sample points by fitting cubic B-spline functions at 
each node, we obtain the smooth power spectrum Psmooth(^)- 

Finally, evaluating the function /bao(^) in each method, we search for the local maxima 
and minima, which we respectively call peaks and troughs to characterize the sound horizon 
scales in BAOs: 

df BAoik) 



= 0. (23) 

peak or trough 



dk 

As a reference, positions of peaks and troughs are computed in each method for the 
linear power spectrum and the resultant numerical values are listed in Table 1. 

4. Baryon Acoustic Oscillations from perturbation theory 

Now we are in a position to discuss the systematic effects on the BAOs by evaluating 
the functions /bao(^)- Consider first the general trends of the systematic influences. In figures 
3, 4 and 5, redshift dependence of the functions /bao is plotted against the wave number for 
three different methods. In each panel, solid lines represent the results from the linear theory 
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Table 1. The linear theory predictions of peak and trough positions in units of /iMpc ^ for our three methods. 



method 


1st peak 


2nd peak 


3rd peak 


4th peak 


1st trough 


2nd trough 


3rd trough 


4th trough 


(i) -P/-Pnw 


0.0669 


0.1221 


0.1792 


0.2377 


0.0415 


0.0944 


0.1500 


0.2075 


(ii) dlnP/dln/c 


0.0533 


0.1074 


0.1638 


0.2216 


0.0823 


0.1368 


0.1933 


0.2514 


(iii) P/Psmooth 


0.0663 


0.1224 


0.1788 


0.2370 


0.0417 


0.0940 


0.1506 


0.2076 



prediction. Note here that the plotted curves from the perturbation results are all restricted to 
the range, k < ki%, where the perturbation theory is safely applied (see Eq. [19]). 

Basically, depending on the characterization methods, the non-linear clustering and the 
redshift-space distortion affect the characteristic scale of BAOs in very different manners. In 
figure 3, deviations from the linear theory prediction become more significant in both real 
and redshift spaces, as the redshift or scale decreases. While the growth of the amplitude 
/bao is clearly seen in real space, the suppression of the amplitude is observed in redshift 
space, together with the overall shifts of acoustic oscillation. The latter effects are simply the 
outcome of the redshift-space distortion. By contrast, systematic effects on the amplitudes 
/bao /bao' shown in figures 4 and 5, seem rather mild. A closer look at the BAO signal 
reveals that the oscillatory features tend to be erased as the redshift and the wavenumber 
increase. This trend has also been seen in A^-body simulations (Seo & Eisenstein 2005; Angulo 
et al. 2007) and modeled by convolving the Gaussian filter with the linear power spectrum 
P^{k) (e.g., Eisenstein et al. 2005). Note that the Gaussian behavior in the disappearance 
of the oscillatory pattern is indeed suggested from the non-perturbative prediction based on 
the renormalized perturbation theory (Crocce & Scoccimarro 2006b). In the present analysis 
using the perturbation theory, the disappearance of the acoustic oscillation mainly comes from 

(22) 

the one-loop terms P^j , showing the monotonic behaviors as function of redshift and wave 
number. 

Turn next to focus on the systematic influences on the location of the peaks and the 
troughs. For this purpose, we define the fractional shift, Ak/k, given by 



k k^ 



(24) 

peak or trough 

and quantify the degree of the positional shifts in each method. Here, the quantities k^"^ and 
f^L ]~ep]~esent the wave number of the peak (or trough) location calculated from the perturbation 
theory and the linear theory, respectively. 

Figures 6, 7 and 8 show the fractional shifts of the peak and the trough positions mea- 
sured from the three different methods, plotted as functions of redshift. In each panel, thick 
and thin lines respectively indicate the fractional shifts for the peaks (from PI to P4) and 
the troughs (from Tl to T4). As a reference, we also plot the observational windows of the 
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planned galaxy redshift survey, WFMOS, depicted as shaded regions. As anticipated from the 
general trends, the systematic changes in the peak and the trough positions become significant 
on small scales and at lower redshift. Also, quite naturally, the magnitude of the shift depends 
on the characterization algorithm. At the redshift around z = 1, the fractional shift of the 
peaks and the troughs reaches or exceeds 3 ~ 4% in cases using /bac while in the method 
/bao' systematic effects are reduced to 1 ~ 2% level. When translating these results to the 
measurement error of the sound horizon scales, the uncertainty in determining the equation of 
state parameter wbe will amount to 12 ~ 16% for the method using /bao' to 4 ~ 8% for 
the method with /bao (^^^ -^^S- !)• the other hand, for the method with /bao) ^ rather 
small value of the fractional shifts (less than 0.5%) is obtained. This is very good news for 
the precision measurement of sound horizon scales. The basic reason for the negligible shifts is 
originated from the definition (22) itself that the non-linear growth and the FOG effect are in- 
corporated into both the numerator and the denominator, by which major sources for the peak 
and the trough shifts are effectively eliminated. As a consequence, the redshift dependence of 
the fractional shifts seen in the real space is almost identical to the one in the redshift space. By 
contrast, the characterization methods with /bao /bao show different redshift dependence 
between real and redshift spaces. Interestingly, the behaviors of the fractional shifts seen in fig- 
ure 6 show some symmetries between peaks and troughs. These systematic trends are related 
to the non-linear corrections arising from the gravitational clustering and the redshift-space 
distortion. We will discussed this issue in some detail in the next section. 

5. A toy model for peak and trough shifts 

The previous section reveals that both the non-linear clustering and the redshift-space 
distortion can lead to the physical shift of the peak and trough positions. However, their 
influences apparently depend on the characterization method for the BAO signals. In this 
section, we discuss how to understand the positional shifts by introducing a simple toy model. 

Our primary goal is a qualitative understanding of the behaviors of the positional shift. 
Let us suppose that the acoustic signature in the linear theory prediction can be described by 
a simple sinusoidal function as 

fBAo{k) = sm{krs + 6o). (25) 

Of course, this is too naive an assumption, but the outcome of the following analysis still keeps 
the essence of our findings. From equation (25), the position of the peaks and the troughs, 
denoted by k^, satisfies the following relation: 

k\, + 6o=[n + ^^7i, (n = 0,l,---). (26) 

Note that the peak (trough) implies that n is even (odd) number. We then add the corrections 
due to the non-linear clustering and the redshift-space distortion, and discuss how the correction 
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induces the positional shifts. Based on figures 3, 4 and 5, the non-hnear corrections may be 
modeled by 

/bao(A;) = e-^^"^" sin(fcr, + 5o) + A{k). (27) 

In the above expression, the Gaussian factor e~''^^^^ represents the smoothing of the acoustic 
oscillations by the non-linear evolution. The function A{k) is assumed to be a monotonic 
function, which mimics the residuals that cannot be absorbed by taking logarithmic derivative 
and/or division by the smooth spectrum. Roughly speaking, monotonically increasing behavior 
of the function A{k) arises from the non-linear growth of gravitational clustering, while the 
monotonically decreasing behavior appears due to the FOG effect in redshift space. Note that 
the parameter A and the function A{k) also depend on the redshift. 

In the present analysis using perturbation theory, the corrections appearing in equation 
(27) should be small and perturbative treatment is always valid. We express the peak and the 
trough position by k^ = k^ + Ak and the shift Ak is treated as small quantity compared to the 
sound horizon scales, l/r^. From the definitions of peaks and troughs, we have (see Eq. [23]), 

dfBAojk) 

peak or trough 

= [-2X^k^sm{k^rs + 6o) + rsCos{k^rs + 6o)}e~'-^''*^^ +A'{k^)=0. (28) 

With a help of the relation (26), the above equation is reduced to the expression for the shift 
Ak, the result of which is summarized as the fractional shift (see Eq. [24]): 

+ : peak 



k k^ \ r„ / k^r 




s 



(29) 

trough 



where we have used the fact that (A/c^)^ <^ k^Vg and k^A' <^ k^Vg. 

From equation (29), systematic changes in the positional shifts may be interpreted as 
a result of the two competing effects. The first term, arising from the Gaussian smoothing 
factor, always makes the position of peaks and troughs move toward smaller k. On the other 
hand, the second term, coming from the monotonic behavior of the residual corrections, affects 
the positional shifts symmetrically: while the peak moves to the high-/c direction, the trough 
is shifted to a small k. The magnitude of these trends will be illuminated more as non-linear 
corrections become important. In figure 9, the role of the two competitive effects are illustrated 
schematically. 

Equation (29) qualitatively explains the behaviors seen in the previous section. In the 
case of the function /bIo' growth or the suppression of amplitudes was significant and the 
smoothing effect of acoustic signature was sub-dominant (Fig. 3). As a result, the peaks and the 
troughs mutually move in an opposite direction, consistent with the toy model (29). A closer 
look at the late-time evolution in redshift space shows somewhat curious behavior that the time 
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evolution of positional shift eventually changes its direction from smaller k to larger k for peaks, 
and from larger k to smaller k for troughs (Fig. 6). Perhaps, this might result from the imbalance 
of the two competitive effects: non-linear growth of gravitational clustering and suppression 
by the FOG effect. Hence, if we allow the sign of A'{k) to change, this is also explained by 
the toy model (29). On the other hand, for the method using /baq ^^"^ /bao' disappearance 
of acoustic oscillation is the major effect (see Figs. 4 and 5). Although the acoustic signature 
seen in the function /baq primarily declined, centered around dlnP/dlnk c:^ —1-7, this does 
not essentially affect the positional shift. As a consequence, the position of the peaks and 
the troughs systematically moves to the low-/c direction. Again, we emphasize the remarkably 
small shift found in the function /bao (FiS- 8)- This implies that the corrections corresponding 
to the A{k) term are completely eliminated. Hence, with the characterization method /baoj 
systematic error in the measurement of sound horizon scale would be greatly reduced, leading 
to an accurate determination of the equation of state parameter wbe- 

6. Conclusions 

We have considered the shifts of the BAO characteristic scales due to the nonlinear 
gravitational and redshift-space distortion effects in a weakly nonlinear regime using one-loop 
correction from perturbation theory. We adopted three different methods to define the BAO 
oscillatory features from the entire power spectrum, and compute the shifts of peak and trough 
locations relative to the purely linear theory predictions in real space. 

In doing so, we presented an analytic toy model to account for the physical reasons for 
the shifts, and showed that one particular method similar to the earlier proposal by Percival 
et al. (2007) is fairly free from the nonlinear and redshift-space distortion effects. In practice, 
the shifts of the first few peak and trough locations defined in the above procedure are at the 
^0.5% level, ensuring ^2% precision in terms of the dark energy parameter wbe, even if one 
uses the linear theory predictions as a standard ruler. Of course the shifts can be accurately 
computed using our methodology as long as the one-loop correction is dominant in the regime 
of interest. 

The result is fairly robust against possible additional effects such as a weakly scale- 
dependent biasing and a running spectral index because they preferentially change the smoothed 
component in the power spectrum that is almost removed from the above procedure. 

Of course the next task is to establish an accurate model that predicts the amplitude of 
the BAO under the nonlinear and redshift-space distortion effects. We suspect that it is very 
challenging given the limitation in both the current theoretical framework and numerical simu- 
lations. In particular, the nonlinear stochastic nature of the galaxy biasing seems problematic 
(e.g., Taruya 2000; Nishimichi et al. 2007). It is still very difficult to incorporate this effect into 
theoretical predictions in any realistic and believable manner. 

In light of this, it may be reasonable at this point to use the BAO scale information 
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exclusively in constraining dark energy, ignoring its amplitude. In practice, constraints on 
cosmological parameters from the use of the BAO scale as a standard ruler hinges on the 
feasibility of the simultaneous fitting of the multiple BAO peaks and troughs. Furthermore, 
two-dimensional (line-of-sight and plane of the sky directions) features in redshift-space may 
improve the accuracy on the BAO scale. We are currently working on a simulation-based study, 
which is necessary to investigate these issues quantitatively. 

We thank M. Takada, A. Nishizawa and E. Reese for useful comments related to the topic 
in the present paper. This work is supported in part by Japan Society for Promotion of Science 
(JSPS) Core-to-Core Program "International Research Network for Dark Energy". T.N, A.S 
and K.Y acknowledge the support from the JSPS Research Fellows. A.T and K.Y are supported 
by a a Grant-in- Aid for Scientific Research from the JSPS (Nos. 18740132, 18540277, 18654047). 

Appendix A. The linear theory prediction of the characteristic scales 

The BAO characteristic scale imprinted in the matter power spectrum is basically the 
sound horizon scale at recombination, rs{z,.ec) (eq.[l]). Depending on the specific definitions 
of peaks and troughs in fc-space that we adopted here, however, their corresponding scales are 
slightly different from the value of equation (1), which has a non-negligible effect in estimating 
the cosmological parameters. The purpose of this Appendix is to clarify the difference in linear 
theory predictions. 

The observed BAO scale in real space is expected to differ slightly from Vslzj-ec) due to the 
residual baryon-photon interaction after recombination {z^^c ~ 1089). Eisenstein & Hu (1998) 
pointed out that a more accurate value is given simply by replacing the z^ec with the drag epoch 



where Rd = R{zd) is the ratio of the baryon to photon momentum densities at = 1019. 
Equation (Al) implies that rs(2rf) = 155Mpc, which is about 5% larger than rs{z^cc) = 148Mpc. 

The characteristic scale in k-space is even more subtle. So let us first model the oscillating 
part of baryon transfer function as 



Eisenstein & Hu (1998) took into account of baryon density perturbation at the drag epoch 
itself, which changes the phase at large scales where the velocity overshoot is not the dominant 
effect. As a result, they found that the phase is approximated as 




(Al) 



Tb{k) oc sin</). 
If we adopt equation (1), the phase is written as 

= krs{Zrcc)- 



(A2) 



(A3) 
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4> = kfs{k), (A4) 

where 

"'^^^ ^ [l + {(3.ojkrs{z,)r]y^' ^^^^ 

/?node = 8.41(fi^/l2)0-435. (A6) 

In these expressions, the characteristic scales in /c— space are defined through 

= ^rn, (A7) 

where m = 5,9, 13, 17, ... for peaks m = 3, 7, 11, 14, ... for troughs from our methods (i) and (iii). 
In contrast, our method (ii) imphes that m = 4,8,12,16,... for peaks m = 6,10,14,18,... for 
troughs. 

Figure 10 compares those theoretical predictions with our peak and trough positions in 
Table 1 (computed using CAMB). The upper panel shows the relation between the phase and 
the wavenumber k; the sohd and dashed lines are calculated using equation (A3) and (A4), 
respectively. The symbols denote our results for peaks and troughs with the three methods 
listed in Table 1. We also plot the fractional deviation with respect to equation (A4) in the 
lower panel. 

The solid and dashed lines suggest that equations (A3) and (A4) are indeed different by 
5% as expected. This is crucial since the difference propagates to ~20% in w^e- The dashed line 
fits our CAMB results (symbols) very well for (p > 57r, which ensures the validity of the formula 
of Eisenstein & Hu (1998). Nevertheless our results for the first few peaks and troughs are 
systematically different from their fitting formula (the lower panel of Figure 10). This would 
simply refiect the difference of our definitions of the peaks and troughs with respect to theirs, 
which one should keep in mind when performing an actual statistical analysis of real datasets. 

Appendix B. Kernels for perturbative solution and one-loop power spectrum 

Here, we briefly summarize the kernels for perturbative solutions presented in equation 
(7) and present the explicit expressions for the one-loop power spectrum (see Eqs.[ll] and [12]). 

In the Einstein-de Sitter universe, the kernel ^i"^(qi, ■ ■ ■ ,qn) satisfles the following re- 
cursion relation (Goroff et al. 1986): 

^i^Hqi) = (i,i), (Bi) 

n-l 

^i'^^ (qi, ■ ■ ■ , qn) = aabin) ^ lUK ki, ka) J'i'") (qi, ■ ■ ■ , q^) 

m=l 

qn— m+1 ; ' ' ' ! ), (B2) 
where k = qi H h q„, ki = qi H h q^, ka = q^+i H h q„, and 

2n + l 2 
3 2n 
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(Jabin) 



{2n + 3){n-l) 



(B3) 



From these relations, one obtains the symmetrized kernels, Fj"^ (density part) and G^""-* (velocity 
part) (e.g., Jain & Bertschinger 1994): 

5 , 1 qi ■q2 /gi , , 2 (qi •q2)^ 



qi,q2j 



7 ' 2 51^2 \q2 ' qij ' 7 qlql 



(B4) 



qi,q2j 



3 ^ 1 qi • q2 /gi ^ ^ 

7 2 gigs U2 ^ g J ^ 7 gfg| 



q2 \ ^4(qi-q2)^ 



(B5) 



Ff(qi,q2,q3) = ^|^^Ff(qi,q2) + 
+cyclic. 



7 k-(qi + q2) ^ 2 Pq3-(qi+q2) 
9 |qi + q2p 9 |qi + q2Pg3 



G'fHqi,q2)| 
(B6) 



(qi,q2,q3) 



l[lk-q3 

613 ql 
+cyclic. 



(qi,q2) + 



1 k ■ (qi + qs) ^ 2 fc^qa ■ (qi + q2) 
3 |qi + q2P 3 |qi + q2Pg| 



(B7) 



The explicit expressions for one-loop correction terms in equations (11) and (12) are 
(e.g., Makino, Sasaki & Suto 1992) 



7^3 /"Oo /* 1 



98(27r)2 

252(27r) 
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(3r + 7x — lOrx 
(l + r2-2rx)^ 



2n2 



-(B8) 



pL(A;) / drP^{kr) 



X 



158 + lOOr^ - 42r^ + ^(r^ - l)^(7r2 + 2) In 



1 +r 



1 — r 



for the power spectrum of density field. 



98(27r)2 Jq 
k^ 



drP^{kr) / dxP^ (kVl + - 2rx 



—r + 7x — 6rx^)^ 
(1 + — 2rx)" 



(B9) 



(,B10) 



^^^^''^^^^ = 84(2^^'^^^/^"^'^^"^ 

X 



12 3 
— - 82 + 4r2 - 6r^ + ^(r^ - lf{r^ + 2) In 



1 +r 



1 — r 



(Bll) 



for the cross power spectrum of density and velocity divergence, and 

98(277)2 

(3r + 7x — lOrx"^) (— r + 7a; — Grx^' 



1/^^ /*00 /* 1 

^^^^''^^^^^ = 98(2^ i) li (Wl + r2-2rx 



X 



(l+r^ — 2rx)2 ' 

;>oo 

-P^(A;) drP^{kr) 

3 

^ - 202 + 56r2 - 30r^ + ^(r^ - l)^(5r2 + 4) In 



(B12) 



252(27r 
24 
2 



1 +r 



1 — r 



(B13) 
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for the power spectrum of velocity divergence. 
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Fig. 3. The power spectrum divided by no-wiggles approximation, /bao(^)' '"^^^ (left) and redshift 
(right) spaces (see Eq.[20]). The solid lines represent the results for the linear power spectrum. The 
others indicate the results for one- loop power spectrum at rcdshifts shown in the panels. The re- 
sults are restricted to the range, k < fci^j, where the perturbation theory is safely applied (Eq.[19]). 
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Fig. 4. Same as figure 3, but for /bao(^) equation (21) 
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Fig. 5. Same as figure 3, but for /baqC'') ™ equation (22) 




Fig. 6. The fractional shifts of peaks (from PI to P4) and troughs (from Tl to T4) of /bao('^)- 
The left (right) panel shows the results in real (redshift) space. The two shaded regions around 
z I and z ~ 3 arc the observational windows of the planned galaxy redshift survey, WFMOS. 
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Fig. 7. Same as figure 6, but for /baqC'') 
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Fig. 8. Same as figure 6, but for /BAo(fc) 




Fig. 9. Schematic figures of the two reasons for the shifts. The upper two curves show /bao(^) in 
case of A'{k) = A = 0, while the lowers represent those for A'{k) > (left) and A > (right). The 
vertical thin dotted lines represent the peak and the trough positions for upper curves, while the short 
vertical solid lines mark those for lower curves. The shifts are given approximately by equation (29). 




Fig. 10. The positions of peaks and troughs predicted from linear theory. 
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